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Abstract: We use dwarf spheroidal galaxies as a tool to attempt to put precise 
lower limits on the mass of the dark matter particle, assuming it is a sterile neutrino. 
We begin by making cored dark halo fits to the line of sight velocity dispersions 
as a function of projected radius (taken from Walker et al. 2007) for six of the 
Milky Way's dwarf spheroidal galaxies. We test Osipkov-Merritt velocity anisotropy 
profiles, but find that no benefit is gained over constant velocity anisotropy. In 
contrast to previous attempts, we do not assume any relation between the stellar 
velocity dispersions and the dark matter ones, but instead we solve directly for the 
sterile neutrino velocity dispersion at all radii by using the equation of state for 
a partially degenerate neutrino gas (which ensures hydrostatic equilibrium of the 
| sterile neutrino halo). This yields a 1:1 relation between the sterile neutrino density 

and velocity dispersion, and therefore gives us an accurate estimate of the Tremaine- 
Gunn limit at all radii. By varying the sterile neutrino particle mass, we locate the 
minimum mass for all six dwarf spheroidals such that the Tremaine-Gunn limit is 
not exceeded at any radius (in particular at the centre). We find sizeable differences 
between the ranges of feasible sterile neutrino particle mass for each dwarf, but 
interestingly there exists a small range 270-280eV which is consistent with all dSphs 
at the l-cr level. 
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1. Introduction 

The lack of success in detecting cold dark matter (CDM) particles |l| emerging 
from supersymmetric extensions to the standard model of particle physics ([[| and 
references therein) has led many authors |, |, |, |, [7|, §, |, 0, |TT], 0, to 
investigate the properties of warm dark matter (WDM) candidates such as sterile 
neutrinos (SNs). WDM might also provide a straight-forward explanation to the 
missing satellites problem [15, [H| and the cusp/core problem [17, 18], both of which 



are of particular significance for the dwarf spheroidals of the Milky Way 

Sterile neutrinos are hypothetical neutral leptons, which don't participate in the 
weak or strong interactions. They can however influence the Universe gravitationally, 
and are also coupled to the active neutrinos through neutrino oscillations. It is 
through these oscillations that the most promising mode of unambiguous detection 
occurs, although no strict limits have been placed on sterile neutrino masses by the 



controversial experiments LSND [|19| and MiniBoone (10 



A thorough study by |2(| attempted to estimate the minimum mass of SNs in 



dwarf spheroidals by taking the average stellar velocity dispersions (VDs), assuming 
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the SNs have the same VD and applying the Tremaine-Gunn inequality pTJ, dis- 



cussed in § |3.3| . In that study, the kinematic data of the dSphs was limited to central 
velocity dispersions (VDs) of a handful of stars, leaving a degeneracy with the ve- 
locity anisotropy. To rectify this general problem, [^2j published line of sight (los) 
velocity dispersions of ~6600 stars for 7 dSphs. This large number of stars (~ 2400 
for Fornax) means we can bin the stars in increments of the projected radius and re- 
duce the degeneracy with the velocity anisotropy. Therefore the dark matter density 
profiles of these dSphs are far better constrained, making studies that predate these 
observations rather obsolete. 

There are a further two flaws in certain previous attempts to set limits on the 
minimum mass of DM particles in dSphs. Firstly, the velocity distribution of the 
DM particles is often erroneously assumed to be the same as that of the stars. 
Realistically, the SNs should exist in a stable equilibrium distribution that satisfies 
hydrostatic equilibrium. In particular, the density of SNs is set by fitting the los VDs 
of the dSphs, but then the equation of state of a partially degenerate neutrino gas 
must be used to compute the unique SN VD for a given SN particle mass, m v . 
Finally, the phase space density found by combining the necessary SN density and 
the VD, must be lower than the TG phase space limit. This can make a considerable 
difference to the deduced minimum mass of SN. 

A further flaw is assuming a priori that the dSph has a cusped dark matter 
halo. As pointed out by |23) and the los VDs cannot distinguish between cored 
or cusped halos yet. There is simply too much freedom in the shape of the halo 
and the velocity anisotropy. The key point here is that cored halos can accomodate 
lower masses of SN than cuspy halos and it is for this reason, rather than seriously 
advocating cored halos, that we use them in our analysis. Cored halos were also 
advocated by p3| and to help explain the large core size needed in the Fornax 
dSph to help explain the survival of the globular clusters in orbit under the influence 
of dynamical friction. 



2. Data and models 



We take the losVDs given in |22| and assume that the stellar samples are not polluted 



by any non-equilibrium stars as suggested in P7| |. In order to minimise the different 
density profiles used in our modelling we apply the same cored density models for the 
DM halos to give new fits to the 3-D light profiles of the dSphs which are excellent 
representations of King models interior to the classical tidal radius (see Fig |I[ Beyond 
the tidal radius, the King models crash towards zero density, but our cored models 
decline with a constant logarithmic slope - which might be more in tune with recent 
observations of the edges of dSphs (e.g. pEfl). The density models are defined as 
follows 
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-7i 

(2-1) 

where the core radius, central density and asymptotic slope are given by r*, p 0ii and 
7i respectively. The subscript "i" corresponds to a "*" for the stellar model or a V 
for the neutrinos. The details for the new stellar density fits are given in Table 1, as 
are the fitted dark halo parameters discussed later. 



Pi(r) = pa, 



1 + 



3. Solving for sterile neutrino equilibrium configurations 
3.1 dSph Jeans analysis 

The density of SNs in dSphs is given by the following steps. Firstly, we must take 
the spherical Jeans equations 

^a 2 (r) + [a(r) + 2/3 (r)] r^a^r) = -g(r) (3.1) 

which solves for the radial velocity dispersion, of(r). For the velocity anisotropy, we 
initially tried an Osipkov-Merritt (OM) model 

0(r) = (3 + 1 (3.2) 

\ r ) 

where (3 , (3^ and rp are free parameters (although -oo < j3 < 1). However, we 
found no benefit from this model over constant anisotropy, although we still solve 
the equations for the OM model. 

The function a(r) is the logarithmic derivative of the 3-D light profile defined in 
Eq U such that 

( \ r dp* -27* 
p m (r) dr 1 + (f) 2 



The density model for the SN halo, given by Eq |2J] also defines the gravity at all 
radii (where we ignore the puny contribution from p*) 

9( r ) = — — / Pv(r)r 2 dr. (3.4) 
r Jo 

From here we use the integrating factor 

J(r) = (r 2 + r 2 )" 7 * (r 2 + r 2 )^r 2/3 °, (3.5) 
which finally allows us to numerically integrate to find the radial velocity dispersions 

POO 

of(r)=r 1 (r) I(r)g{f)dr. (3.6) 
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We then cast the radial velocity dispersion into the line of sight to compare with 
the observations, and minimise the x 2 - The best fit SN density parameters, p 0>v , r v 
and velocity anisotropy are given in Table 1, where in order to minimise the freedom 
of our models, we force 7^ = 1.5. The error bars on the 3 free parameters correspond 
to the differences the parameters would need to have from the best fit in order to 
raise the x 2 by 3.5. As expected, this is a large parameter space, because changes 
in the combination of p 0) „ and r v can be counter balanced by the variable velocity 
anisotropy. Fortunately, the shape of the losVD profile keeps these combinations 
finite. Furthermore, although the individual ranges in p ov and r v are vast, combined 
measures like the enclosed masses (for instance M 30 o in Table 1) are much less ex- 
tensive. In the middle panel of Fig 1, the best fit model is given by the solid line 
and the dashed lines are the 1-cr errors on the velocity anisotropy. 

3.2 Sterile neutrino hydrostatic equilibrium 

Now we have a density of SNs in each dSph fixed by the stellar losVDs: at least it 
is the lowest density exactly required for hydrostatic equilibrium of the stars in the 
dSph, that has centrally isotropic orbits and fits the data well. In addition to this 
constraint, there is an equation for hydrostatic equilibrium of the SNs themselves. 
This invokes the equation of state of neutrinos and gives us a relation that defines 
the density and pressure of a partially degenerate neutrino gas, coupled via the 
hydrostatic equilibrium relation 4pP v = —PvQij')- 

The exact procedure to find the velocity dispersion required for hydrostatic equi- 
librium of the SN halo is rather lengthy, but is given in its full glory in p9[ . The 
end result is that there is a one to one relation between the SN density, p u , and 
gravity (which is defined by p v in Eq |3~5|) and the SN velocity dispersion, with the 
SN particle mass the only tunable parameter. 

3.3 Phase space constraints 

We now have a unique correlation between the density of SNs and their velocity 
dispersions. This is important because these two variables define the phase space 
distribution of the SNs, to which there is a fundamental limit. 

Liouville's theorem states that (in the absence of encounters) flow in phase space 
is incompressible and that each element of phase-space density is conserved along 
the flow lines. However, this only applies to the fine-grained phase space density 
of an infinitesimal region. Rather, for the observable, which is the coarse-grained 
(macroscopic) phase-space density, we simply must not exceed the maximum of the 
fine-grained one. 

Thus, the SN phase-space density must not increase during collapse, from its 
starting value of ^g v h~^m'l (which is half the Pauli limit), to its current maximum 
value of py(r)[27T(Ty(r)] _1 - 5 (where we assume the velocity distribution is a Gaus- 
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sian with dispersion a u ). This limit is called the Tremaine-Gunn (TG) limit and 
rearranging it in terms of the critical density for a given SN particle mass gives us 



in,, \ I cr„(r) 1 ''' 



The important thing to bear in mind here is that this is a more accurate value of the 
TG limit because it is calculated from the unique, derived SN VD exactly necessary 
for hydrostatic equilibrium and not by assuming some relation between the stellar 



VD and the sterile neutrino VD, as is often the case [30, 31, Ell. The only assumption 



is that the velocity distributions are locally Maxwellian everywhere. 

This defines the TG limit as a function of radius, since we know the exact value 
of the VD at all radii from solving for the SN hydrostatic equilibrium. Our final step 
is to vary the mass of the SN such that the TG limit is reached exactly at the centre. 
For the best fit parameters, the minimum SN mass is given in Table 1 and the errors 
are propogated by finding the minimum SN mass needed when the extreme values 
of p 0jU and r u are used (j3 plays no part). 



4. Results 

In Fig |l| we present the densities, VDs and enclosed mass profiles for the dark halos 
and stars for the six dSphs. In the left hand panels we plot the SN density (solid 
line type) and the TG limit of the dSph (dashed) against radius. In addition to this, 
we show the King model densities (with unity M/L) along with the Eq 1 density 
models (with parameters given in Table 1). In the middle panels we plot both the 
derived SN (dotted) VD and stellar losVD (solid) which is fitted to the data points 
taken from p2] . In the right hand panels, we plot each dSph's enclosed SN dark halo 
(solid) and stellar (dashed) masses. 

It is encouraging that the velocity dispersions of the SNs are similar to the 
losVDs of the stars, and that they are smooth. It is worth mentioning here that 
the equilibrium SN velocity dispersion (and hence TG limit) depends critically on 
the slope the SN density profile. As can be seen in the left hand panels of Fig |T], 
the TG limit remains fairly constant at all radii. Our SN halos, although cored in 
name, continue to increase slowly towards the centre in all cases. It is this behaviour 
that forces the TG limit to fall slightly (a factor of 2 or so only in these cases) from 
slightly outside the centre. If the density increases relatively rapidly, the phase space 
density increases towards the TG limit (also towards the Fermi limit). As it does so, 
the classical pressure required for hydrostatic equilibrium (P v , c oc p v ^1) is gradually 
replaced by Fermi pressure (P u j oc pi! 3 ) causing the VD to decrease, and hence the 
TG limit to decrease. Since we cannot decrease the SN density any further (without 
sacrificing our fits to the stellar losVDs), if the SN particle mass is too small, then 
the TG limit will dip below the required density and will be unfeasible. 
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Table 1: Here we list all the parameters and results pertaining to the fits of our cored sterile neutrino halos to the six dSphs 
and the best fit light profiles for each dSph. Also given are the sterile neutrino mass within 300 pc and the minimum mass of sterile 
neutrino particle that we calculate can exist in hydrostatic equilibrium. These are compared to the traditional Trcmainc-Gunn limits 
(Madscn 1991). In each case, the velocity anisotropics arc constant. The distances and luminosities of the dSphs arc all taken from 
Table 1 of Angus (2008), where the original references are given. Error bars on the three free fit parameters are at the 1-cr error and 
these propogatc simply into the minimum neutrino mass and -M3Q0' 
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This can be thought of another way: the more cuspy the SN density needs to be, 
the larger the SN particle mass must be to evade the TG limit. It is by increasing 
the SN particle mass, that we increase the TG limit at all radii and choose the mass 
that causes the TG limit to drop to the required density only at the centre. 

The minimum SN masses required for our six dSphs are given in Table 1 and 
are plotted against dSph luminosity in Fig 2. Interestingly, although the two Leo 
dSphs provide the extremes of the preferred particle mass (160^^5 and 460+24° f° r 
Leo I and II respectively), the other four lie in a relatively tight range between 250 
and 380 eV. Surprisingly, there exists a tiny 10 eV interval (270-280 eV) which all 
six dSphs are consistent with at 1-a. This minimum SN particle mass is significantly 
smaller, in all but the case of Sextans, than the traditional estimates of the TG limit 



(given in Table 1) found using Eq 2 of |32| . 

Our approach also improves the limits on the lowest allowed SN mass because 
we do not need to assume a relation between the stellar losVD and the velocity 
dispersion of the SNs (we solve directly for it). In five of the six cases, the SN VD is 
20-100% larger than the stellar losVD, leaving the unusual case of Leo I which prefers 
a very large core - meaning a lower mass SN with a larger velocity dispersion can be 
accomodated. Furthermore, since we use cored neutrino halo models, we allow less 
massive neutrinos than cuspier halos would allow. 

[ |33|1 proposed the common mass scale where each dSph dark halo should have 
~ 10 7 M Q within 300 pc, the ultra-faint dSphs |34|, |35|, [36|, [37], ^ discovered by the 
SDSS should pose a more stringent test for us. However, their small samples of stars 
mean we cannot reliably use our method to determine an equilibrium configuration 
for the SNs. In addition, questions have been raised by |39| and p7j about how well 
interlopers have been dealt with in these systems. Therefore, our limits should be 
the strongest and most reliable ones possible. 

On the same topic, our cored halo models offer a sanity check on the cuspier halo 
models used to constrain the mass within 300 pc (M 30 o) for all dSphs by fl3~3f . Like 
j33fl , our M 30 os do not appear to be correlated with satellite luminosity, although the 



preferred common mass scale is very slightly (10-20%) higher. 
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5. Summary and Conclusion 



We have taken the data of on the line of sight velocity dispersions of six dwarf 
spheroidal galaxies (dSphs) and fitted cored sterile neutrino halos assuming some 
anisotropic orbits for the stars. From this we solved for the velocity dispersion of 
the dark halos at all radii by solving the equation of state for a partially degenerate 
neutrino gas. This allowed us to infer the lower limit on the dark matter particle 
mass (assuming it is a sterile neutrino) and found it is in the range 160-460 eV. This 
limit is lower than estimates achieved with the traditional Tremaine-Gunn approach 
and intriguingly, there exists a range of minimum sterile neutrino mass (270-280 eV) 
which is consistent with all six dSphs at the 1-a level. 
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Figure 1: Each row refers to a specific dwarf spheroidal (identified in the left hand panel), 
whereas the columns are properties of that dwarf spheroidal. In the left hand panels we 
plot the 3-D density against radius for the best fit sterile neutrino halo (solid), and the 
corresponding TG limit (dashed). In addition we plot the King model stellar density 
(dotted, with M/L=l) and overplot the best fit cored density model of Eq 1. In the 
central panels we plot the inferred 3-D velocity dispersion of the sterile neutrinos (dashed) 
required for their hydrostatic equilibrium, and the projected stellar line of sight velocity 
dispersions (solid) - which are fitted to the data points (asterix symbols) from Walker et 
al. (2007). The dotted lines are the upper and lower 1-a errors given by the extremes of 
velocity anisotropy. In the right hand panels we plot enclosed sterile neutrino halo (solid) 
and stellar (dotted) masses against radius. 
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Figure 2: The left hand figure shows the minumum sterile neutrino masses plotted against 
their luminosities. The thin dashed strip shows the range of neutrino masses that are 
consistent with all six dSphs at the 1-a level. In the right hand figure, we plot M300 
against luminosity. Our six dSphs seem to prefer a slightly higher value of M300 than 
the value deduced in Strigari et al. (2008), perhaps only due to the different dark matter 
density profiles and velocity anisotropies. 



